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Abstract 

Rheumatoid arthritis (RA) is a chronic systemic inflammatory disease. Certain HLA-DRBl "shared-epitope" alleles are 
reported to be positively associated with increased RA susceptibility, whereas some of the other alleles may be negatively 
associated. However, studies on the latter are rare. Here, we focus on the protective effects of DRBl alleles in Japanese RA 
patients in an association study. Relative predispositional effects (RPE) were analyzed by sequential elimination of carriers of 
each allele with the strongest association. The protective effects of DRBl alleles were investigated in patients stratified 
according to whether they possessed anti-citrullinated peptide antibodies (ACPA). The DRB1*13:02 allele was found to be 
negatively associated with RA {P = 4.59x10"^° corrected P (Pc) = 1.42x1 0"^ odds ratio [OR] 0.42, 95% CI 0.32-0.55, P 
[RPE] = 1.27x10"*); the genotypes DRBl *04:05/*l 3:02 and *09:01/*13:02 were also negatively associated with RA. The 
protective effect of *73:02 was also present in ACPA-positive patients (P = 3.95x10"**, Pc = 1.22xl0"^ OR 0.42, 95%CI 0.31- 
0.58) whereas *15:02 was negatively associated only with ACPA-negative RA (P = 8.87xl0"^ Pc = 0.0026, OR 0.26, 95%CI 
0.12-0.56). Thus, this study identified a negative association of DRB1*13:02 with Japanese RA; our findings support the 
protective role of DRB1*13:02 in the pathogenesis of ACPA-positive RA. 
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Introduction 

Rheumatoid arthritis (RA) is a chronic systemic inflammatory 
disease that affects about 1 % of the population. Its pathogenesis is 
multifactorial and disease susceptibility is associated with genetic 



and environmental factors [1,2,3]. Human Leukocyte Antigen 
(HLA) alleles are associated with RA in most ethnic groups and 
represent the strongest genetic risk factors for the disease. Most 
reports are of HLA-DRBl alleles positively associated with RA 
susceptibility. A conserved amino acid sequence at position 70-74 
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(QKRAA, RRRAA, or QRRAA) in the HLA-DRP chain is 
shared between the RA susceptibility-associated DRBl alleles; this 
was designated the "shared epitope" (SE) [4]. A gene dosage effect 
was noted in the associations of HLA-DRBl alleles with 
susceptibility to RA in that homozygosity for susceptibility alleles 
does confer higher disease risk than heterozygosity for these alleles. 

The presence of anti-citruUinated peptide antibodies (ACPA) is 
associated with RA with higher specificity than rheumatoid factor; 
thus, ACPA is thought to play some role in the pathogenesis of 
RA, especially as SE alleles are strongly associated with ACPA- 
positive RA but only relatively weakly with ACPA-negative RA 
[5]. Several studies have found that DRB1*04:01 and *04:05, both 
SE alleles, were mainly associated with RA in European and East 
Asian populations, respectively. 

As well as associations with disease susceptibility, some DRBl 
alleles are reported to be negatively associated with RA. An amino 
acid sequence (DERAA) at position 70-74 [6], isoleucine at 
position 67 (167) [7], aspartic acid at position 70 (D70) [8], or a 
conserved amino acid sequence at position 71—74 (SI; ARAA or 
ERAA) [9,10] in the HLA-DRfi chain seem to be protective in 
European populations. It was also reported that DRB1*13 alleles 
are negatively associated with ACPA-positive and -negative RA in 
European populations [11]. A meta-analysis concluded that 
DRB1*13:01 was protective against ACPA-positive RA in 
European populations [12]. However, there are very few studies 
on the protective effects of DRBl alleles in Japanese patients, 
although reduced frequencies of some DRBl alleles have been 
reported in Asian RA [13,14,15,16,17]. In this study, we focus on 
the protective effects of HLA-DRBl alleles in Japanese RA patients 
with or without ACPA. 

Materials and Methods 

Patients and controls 

RA patients (n = 1480) were recruited at Sagamihara Hospital, 
Tama Medical Center, Nagoya Medical Center, Nagasaki 
Medical Center, Yokohama Minami Kyosai Hospital, Kumamoto 
Center for Arthritis and Rheumatology, Miyakonojo Hospital, 
Niigata Rheumatic Center, and Hyogo College of Medicine. Of 
these 1480 RA patients, 919 were ACPA-positive and 110 were 
ACPA-negative. ACPA data were not available for the remaining 
45 1 patients. Healthy controls (n = 800; mean age ± SD, 
36. 7± 10.7 years, 238 male [30.1%]) were recruited at Sagamihara 
Hospital and University of Tokyo, or by the Pharma SNP 
Consortium (Tokyo, Japan) [18]. All patients and healthy 
individuals were native Japanese living in Japan. All patients with 
RA fulfilled the 1987 American College of Rheumatology criteria 



for RA [19]. Rheumatoid factor and ACPA were detected using 
the N-latex RF kit (Siemens Healthcare Diagnostics, Miinchen, 
Germany) and the Mesacup-2 test CCP (Medical & Biological 
Laboratories, Nagoya, Japan), respectively. This study was 
reviewed and approved by the Research Ethics Committees of 
each participating institute: Nagasaki Medical Center Research 
Ethics Committee, Yokohama Minami Kyosai Hospital Research 
Ethics Committee, Tama Medical Center Research Ethics 
Committee, University of Tsukuba Research Ethics Committee, 
Miyakonojo Hospital Research Ethics Committee, Kumamoto 
Center for Arthritis and Rheumatology Research Ethics Commit- 
tee, Niigata Rheumatic Center Research Ethics Committee, 
Hyogo College of Medicine Research Ethics Committee, and 
the University of Tokyo Research Ethics Committee. Written 
informed consent was obtained from aU study participants. This 
study was conducted in accordance with the principles expressed 
in the Declaration of Helsinki. 

Genotyping 

Genotyping of HLA-DRBl was performed by a polymerase 
chain reaction technique using sequence-specific oligonucleotide 
probes (WAKFlow HLA typing kits, Wakunaga, Hiroshima, 
Japan), using a Bio-Plex 200 system (Bio-Rad, Hercules, CA), or 
using MPH-2 High Resolution HLA typing kits (Wakunaga) for 
four-digit allele typing. The following DRBl alleles contain the SE 
[4]: *01:01, *04:01, *04:04, *04:05, *04:10, *10:01, *14:02, and 
*14:06. DRBl aUele groups, D70, 167, SI, and DERAA, were 
reported to be protective in European populations [6,7,8,9,10]; the 
protective effects of these allele groups in Japanese were validated 
in this study. DRBl aUeles containing D70 [8] are *07:01, *08:02, 
*08:03, W8:()9, *08:23, *11:01, *11:06, *12:01, *12:02, *12:05, 
*13:01, *13:02, "13:07, * 14:03, * 14:1 2, and * 16:02. DRBl aUeles 
containing 167 [7] are *07:01, *08:03, *08:23, *12:01, *12:05, 
*13:01, *13:02, *14:45, *15:01, *15:02, and "'15:11. DRBl alleles 
containing DERAA [6] are the same as DRB1*13 (i.e. "13:01, and 
*13:02). FinaUy, DRBl aUeles containing SI [20] are "13:01, 
*13:02, *15:01, and *15:02. Results oi DRBl genotyping for some 
of the healthy controls were reported previously [14]. Some of the 
RA patients were also included in another study which reported 
on susceptibility effects for interstitial lung disease or positivity for 
autoantibodies [21,22,23]. HLA-DRBl genotype of each subject 
was not deposited in publicly available resources. 

Statistical analysis 

The exact tests for deviation from Hardy- W einberg equilibrium 
were conducted by the Markov chain method under the condition 



Table 1. Characteristics of the RA patients studied. 







RA 


ACPA(+) RA 


ACPA(-) RA 


P 


Number 


1480 


919 


110 




Mean age, years (SD) 


63.9 (12.2) 


63.7 (12.2) 


63.4 (12.3) 


0.8582* 


Male, n (%) 


272 (19.0) 


171 (18.7) 


21 (19.3) 


0.8969 


Age at onset, years (SD) 


49.3 (14.4) 


49.2 (14.2) 


50.1 (16.7) 


0.6092* 


Steinbrocker stage III and IV, n (%) 


560 (37.8) 


521 (56.7) 


45 (40.9) 


0.0703 


Rheumatoid factor positive, n (%) 


1002 (67.7) 


826 (89.9) 


40 (36.4) 


9.39x10"" 



RA: rheumatoid arthritis, ACPA: anti-citrullinated peptide antibody, ACPA(+): ACPA-positive, ACPA(— ): ACPA-negative. Association was tested by Fisher's exact test using 
2x2 contingency tables or Student's t-test. *Student's t-test was employed. 
doi:l 0.1 371/journa!.pone.0099453.t001 
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of 10000 each of dememorization, batches, and iterations per 
batch (Genepop on the web; http://genepop.curtin.edu.au/) [24]. 
DrfiFerences of allele carrier frequencies, genotype frequencies or 
amino acid residue carrier frequencies were analyzed by Fisher's 
exact test using 2x2 contingency tables. In order to estimate the 
protective effects of alleles in multi-aUelic locus on individuals for 
RA, differences of allele carrier frequencies, or amino acid residue 
carrier frequencies were analyzed under the dominant model. 
Adjustment for multiple comparisons was performed using the 
Bonferroni method. Pc values were calculated by multiplying the P 
value by the number of alleles or amino acid residues tested. 

Alleles with low carrier frequencies in RA patients may not be 
detectably protective because predisposing SE alleles with higher 
carrier frequencies could obscure their influence. To investigate 
the protective effects of HLA alleles, relative predispositional 
effects (RPE) were analyzed by sequential elimination of carriers of 
each allele with the strongest association [25]. In order to obtain 
an accurate estimate of the effects of alleles other than SE, analyses 
of these alleles in RA patients were also stratified in the following 
manner: For SE-negative subjects, the effect in ^'A/A^' and "A/ 
other than SE or A" genotype groups was investigated using 
"other than SE or yl/other than SE or A" genotype group as the 
reference. For SE-positive subjects, the effect of "SE/^" genotype 
group was analyzed using "SE/other than A" genotype group as 
the reference. The protective effects of the *13:02 allele were 
confirmed in the presence of predisposing allele "5". The effect in 
"B/ *13:02" genotype group was investigated using the "iJ/ other 
than *13:02" genotype group as the reference. The protective 
effects of the *15:02 allele were confirmed in the analysis of "B/ 
*15:02'' using the "i?/ other than *15:()2" genotype group as the 
reference in the same manner. 

Results 

Characteristics of RA patients 

Characteristics of ACPA-positive [ACPA(+)] and ACPA-nega- 
tive [ACPA(— )] RA patients are given in Table 1. The proportion 
of rheumatoid factor-positive patients in the ACPA(+) group was 
higher than in ACPA(— ) RA. There were no significant differences 
in terms of mean age, percentage of males, age at onset, or 
Steinbrocker stage [26] between ACPA(+) and ACPA(-) patients. 



Reduced HLA-DRB1*13:02 allele carrier frequency in 
Japanese RA 

n = HLA-DRBl genotyping was performed in 1480 RA patients and 

^ S 800 healthy controls to compare HLA allele carrier frequencies 

S '-^ (Table 2). No deviation from Hardy-Weinberg equilibrium was 

I .2 observed in the controls (/"= 0.6329), though a deviation was 

m S detected in the RA patients (P<0.0001). A strong positive 

^. -3 association between the frequency of DRB1*04 and RA 

I I (P= 1.00x10"^^ Corrected P [Pc] = 1.31 xlO"^', odds ratio 

[OR] 2.40, 95% confidence interval [CF] 2.01-2.86, Table 2) 
was confirmed. Additionally, DRB1*13 (i.e. the DERAA allele 
group) was found to be negatively associated with RA 
(P=4.69xl0"", Pc = 6.10xl0"'", OR 0.41, 95% CI 0.31- 
0.53). The D70, 167, and SI allele groups were also negatively 
associated with RA (D70: P= 1.15x10"^", OR 0.43, 95% CI 
0.36-0.52; 167: P=3.67xl0"'^ OR 0.52, 95% CI 0.44-0.62; SI: 
P=2.51xl0"^ OR 0.59, 95% CI 0.49-0.70). Finally, a 
predisposing association was confirmed between SE and RA 
(P=2.35xlO"*\ OR 3.58, 95% CI 2.99-4.29). 

We further explored associations between these DRBl alleles 
and RA by high-resolution typing, using RPE testing [25] 
(Table 2). RPE were analyzed by sequential elimination of carriers 
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Table 3. HLA-DRBl a 


lele carrier frequency in 


RA patients and controls in su 


bjects stratified for the presence of SE. 








RA (n= 1480) 


Control (n = 


= 800) 


P 


OR 


95%CI 


'03 


SE negative 


1 (0.2) 


0 (0.0) 




0.4785 


3.28 


(0.13-80.65) 




SE positive 


1 (0.1) 


0 (0.0) 




1 .0000 


0.91 


(0.04-22.52) 


•04 other than SE 


SE negative 


46 (10.3) 


81 (16.7) 




0.0055 


0.58 


(0.39-0.85) 




SE positive 


54 (5.2) 


34 (10.8) 




0.0010 


0.45 


(0.29-0.71) 


•07 


SE negative 


5 (1.1) 


5 (1.0) 




1 .0000 


1.09 


(0.31-3.79) 




SE positive 


5 (0.5) 


2 (0.6) 




0.6680 


0.76 


(0.15-3.93) 


•08 


SE negative 


114 (25.6) 


140 (28.9) 




0.2703 


0.85 


(0.64-1.13) 




SE positive 


74 (71) 


41 (13.0) 




0.0018 


0.51 


(0.34-0.77) 


•09 


SE negative 


230 (51.7) 


161 (33.2) 




1.44x10"" 


2.15 


(1.65-2.81) 




SE positive 


193 (18.6) 


52 (16.5) 




0.4051 


1.16 


(0.83-1.62) 


•77 


SE negative 


17 (3.8) 


29 (6.0) 




0.1336 


0.62 


(0.34-1.15) 




SE positive 


23 (2.2) 


4 (1.3) 




0.3635 


1.77 


(0.61-5.15) 


•72 


SE negative 


63 (14.2) 


66 (13.6) 




0.8496 


1.05 


(0.72-1 .52) 




SE positive 


80 (7.7) 


21 (6.7) 




0.6248 


1.17 


(0.71-1.93) 


•7J 


SE negative 


48 (10.8) 


102 (21.0) 




2.43x10"^ 


0.45 


(0.31-0.66) 




SE positive 


64 (6.2) 


32 (10.2) 




0.0235 


0.58 


(0.37-0.91) 


*;4 other than SE 


SE negative 


63 (14.2) 


89 (18.4) 




0.0918 


0.73 


(0.52-1 .04) 




SE positive 


72 (7.0) 


33 (10.5) 




0.0537 


0.64 


(0.41-0.98) 


•75 


SE negative 


182 (40.9) 


207 (42.7) 




0.5950 


0.93 


(0.72-1.21) 




SE positive 


223 (21.5) 


55 (17.5) 




0.1305 


1.30 


(0.94-1.80) 


•76 


SE negative 


10 (2.2) 


9 (1.9) 




0.8174 


1.22 


(0.49-3.02) 




SE positive 


7 (0.7) 


6 (1.9) 




0.0900 


0.35 


(0.12-1.05) 


* 13:01 


SE negative 


1 (0.2) 


7 (1.4) 




0.0712 


0.15 


(0.02-1.26) 




SE positive 


4 (0.4) 


1 (0.3) 




1 .0000 


1.22 


(0.14-10.94) 


* 13:02 


SE negative 


47 (10.6) 


95 (19.6) 




0.0001 


0.48 


(0.33-0.71) 




SE positive 


60 (5.8) 


31 (9.8) 




0.0148 


0.56 


(0.36-0.89) 


D70 


SE negative 


242 (54.4) 


320 (66.0) 




0.0004 


0.61 


(0.47-0.80) 




SE positive 


261 (25.2) 


114 (36.2) 




0.0002 


0.59 


(0.45-0.78) 


167 


SE negative 


301 (67.6) 


371 (76.5) 




0.0027 


0.64 


(0.48-0.86) 




SE positive 


390 (37.7) 


130 (41.3) 




0.2614 


0.86 


(0.67-1.11) 


SI 


SE negative 


217 (48.8) 


288 (59.4) 




0.001 2 


0.65 


(0.50-0.84) 




SE positive 


287 (27.7) 


87 (27.6) 




1.0000 


1.01 


(0.76-1.33) 


D70 other than *13:02 


SE negative 


212 (47.6) 


253 (52.2) 




0.1892 


0.83 


(0.64-1 .08) 




SE positive 


201 (19.4) 


83 (26.3) 




0.0092 


0.67 


(0.50-0.90) 


167 other than *13:02 


SE negative 


280 (62.9) 


310 (63.9) 




0.7852 


0.96 


(0.73-1.25) 




SE positive 


330 (31.9) 


99 (31.4) 




0.8904 


1.02 


(0.78-1.34) 


SI other than '13:02 


SE negative 


182 (40.9) 


213 (43.9) 




0.3536 


0.88 


(0.68-1.15) 




SE positive 


227 (21.9) 


56 (17.8) 




0.1147 


1.30 


(0.94-1 .80) 



RA: rheumatoid arthritis, SE: Shared epitope, OR: odds ratio, CI: confidence interval. Allele carrier frequencies are shown in parenthesis {%). Association was tested by 
Fisher's exact test using 2x2 contingency tables. SE negative: "A/A" or "yA/other than SE or A" vs. "other than SE or yA/other than SE or A". SE positive: "SE/>A" vs. "SE/ 
other than A". Allele groups SE, D70, 167, and SI were as defined in the Materials and Methods section. 
doi:1 0.1 371 /journal.pone.0099453.t003 



of each allele with the strongest association (Table 2, right 
column). The prime strongest association was between 
DRB1*04:05 and RA (P= 1.41 xlO"'^^ Pc = 4.37 x 10"'", OR 
3.31, 95% CI 2.73-4.01). Thus, a second round of comparisons 
was conducted after the elimination of DRB1*04:05 carriers, 
revealing the next strongest association to be between DRBl* 13:02 
and RA {P= 1.27 x 10"*^, Pc = 3.68x lO""^). A third round after the 
elimination of both DRB 1*04:05 or * 13:02 carriers now showed 



the strongest association of RA with DRB1*04:01 (P= 0.0002, 
Pc = 0.0065). Further rounds after elimination oi DRB1*04:05, 
^13:02 and *04:01 carriers revealed associations between the 
remaining DRBl alleles and RA, particularly for DRB1*09:01 
(P=7.32xlO"\ ft = 0.0020), *01:01 (P= 3.75xlO"\ 
Pc = 0.0010), *14:06 (P= 0.0041, Pc = 0.0995), *10:01 
(P= 0.0128, Pc = 0.2936), *04:10 (P= 0.0109, Pc = 0.2399), and 
*11:01 (P= 0.0236, ft = 0.4948). The results from association 



PLOS ONE I www.plosone.org 



5 



June 2014 | Volume 9 | Issue 6 | e99453 



DRB1*13:02 Protects against RA 



Table 4. HLA-DRBl genotype frequency in RA patients and controls. 







RA (n = 1480) 


Control (n = 800) 


P 


OR 


95%CI 


*04:05/*13:02 


45 (6.1) 


21 (11.4) 


0.0168 


0.51 


(0.29-0.87) 


»04:01/*13:02 


6 (7.1) 


2 (11.8) 


0.6190 


0.58 


(0.11-3.14) 


*09:01/*13:02 


12 (2.8) 


21 (9.9) 


0.0004 


0.27 


(0.13-0.55) 


'01:01 /'1 3:02 


5 (2.4) 


3 (3.6) 


0.6917 


0.65 


(0.15-2.79) 



RA: rheumatoid arthritis, OR: odds ratio, CI: confidence interval. Allele carrier frequencies are shown in parenthesis {%). Association was tested by Fisher's exact test 
using 2x2 contingency tables. Comparison: "B/*13:02" vs. 6/other than *13:02". 
doi:1 0.1 371 /journal.pone.0099453.t004 



Studies under the recessive and the allele models were represented 
in Table SI and S2, respectively. Similar tendencies were observed 
in these analyses. We therefore focused on the DRBl allele with 
the most signilicandy reduced allele carrier frequency, namely 
DRBl *1 3:02. 

Protective effects of the *73;02 allele against RA in both 
SE-positive and -negative subjects 

In order to obtain an accurate estimate of the effects of alleles 
other than SE, associations were estimated in subjects stratified 
into those with or without SE (Table 3). Although DRB1W9 
(P= 1.44x10"^ OR 2.15, 95% CI 1.65-2.81) predisposes to RA 
in SE-negative people, *04 other dian SE {*04:03, *04:06, *04:07: 
SE negative, P= 0.0055, OR 0.58, 95% CI 0.39-0.85; SE positive, 
P= 0.0010, OR 0.45, 95% CI 0.29-0.71), *13{*13:01, *13:02: SE 
negative, P=2.43xl0"\ OR 0.45, 95% CI 0.31-0.66; SE 
positive, P= 0.0235, OR 0.58, 95% CI 0.37-0.91), and D70 (SE 
negative, P= 0.0004, OR 0.61, 95% CI 0.47-0.80; SE positive, 
P= 0.0002, OR 0.59, 95% CI 0.45-0.78) were negatively 
associated with RA in both SE-positive and -negative individuals. 
DRB1*08 (P= 0.0018, OR 0.51, 95% CI 0.34-0.77) aUeles were 
negatively associated with RA in SE-positive people. 167 
(P= 0.0027, OR 0.64, 95% CI 0.48-0.86) and SI (P= 0.0012, 
OR 0.65, 95% CI 0.50-0.84) alleles were negatively associated 
with RA in SE-negative subjects. However, D70 alleles other than 
*13:02 were negatively associated with RA in SE-positive 
(P= 0.0092, OR 0.67, 95% CI 0.50-0.90) but not in SE-negative 
individuals. 167 aUeles other than *13:02 and SI aUeles other than 
*13:02 did not have any negative associations. These data suggest 
that the negative associations of D70, 167 and SI aUeles with RA 
in SE-negative subjects were mainly mediated by *13:02, although 
the negative association of D70 in SE-positive people was due to 
*08 alleles. Thus, * 13:02 was negatively associated with RA in SE- 
negative people and relatively weakly also in SE-positive subjects. 

The protective effects of the *13:02 allele were analyzed in the 
presence of predisposing alleles (Table 4). Although *04:05 and 
*09:01 are positively associated with RA in Japanese, the risk of 
disease in people carrying these aUeles was decreased in 
heterozygotes also carrymg *13:02 [*04:05: P= 0.0168, 
OR = 0.51, 95% CI 0.29-0.87; *09:01: P=0.0004, OR = 0.27, 
95% CI 0.13-0.55). 

Although the age at RA onset in *04:05 aUele carriers was lower 
than in non-carriers (mean age ± standard deviation [SD] [years] , 
carriers vs. non-carriers, 48.2± 13.8 vs. 50.5± 14.9, P= 0.0070) the 
age at onset in people with *13:02 or *01:01 was higher than in 
non-carriers (53.8±14.0 vs. 48.9±14.4, P= 0.0027, and 
52.9±13.3 vs. 48. 7± 14.5, P= 0.0021, respectively) (Table S3). 



Protective effects of *13:02 against ACPA(-i-) RA and 
*/5;02 against ACPA(-) RA 

Predisposing effects of the *04:05 aUele were confrrmed in 
ACPA(+) RA (Table 5, P= 3.64xl0"'^\ Pc = 1.13 x 10"'" OR 
3.59, 95% CI 2.91-4.42), whereas DRB1*13:02 was negatively 
associated with ACPA(+) RA (P=3.95xlO"^ Pc=1.22xl0"^ 
OR 0.42, 95% CI 0.31-0.58). The DERAA aUele group was stiU 
negatively associated with RA even when only ACPA(+) patients 
were considered (P=2.05xlO"^ OR 0.40, 95% CI 0.29-0.54). 
D70, 167, and SI were also negatively associated with ACPA(+) 
RA (D70: P= 5.78 x 10"^', OR 0.39, 95% CI 0.32-0.48; 167: 
P=3.66xl0"'^ OR 0.50, 95% CI 0.42-0.61; SI: P=2.31xl0"^ 
OR 0.57, 95% CI 0.47-0.70), and the predisposing association 
was confirmed between SE and ACPA(+) RA (P= 1.16xl0"*^ 
OR 4.41, 95% CI 3.59-5.41). 

The predisposing association was also confirmed between SE 
and ACPA(-) RA (P= 0.0229, OR 1.60, 95% CI 1.07-2.38), 
albeit weakly. A tendency towards a positive association of *04:05 
and *14:54 with ACPA(-) RA was observed {*04:05: P= 0.0126, 
Pc = 0.3667, OR 1.75, 95% CI 1.1 5-2.69; *i^:,5^: P= 0.0202, OR 
2.25, Pc = 0.5861, 95% CI 1.17-4.32). On the other hand, 
DRBl *1 5:02 was negatively associated with ACPA(-) RA 
(P=8.87xl0"\ Pc = 0.0026, OR 0.26, 95% CI 0.12-0.56). 

We next examined associations of alleles other than SE with 
ACPA(+) and ACPA(— ) RA stratified by the presence or absence 
of SE (Table 6). This analysis showed that DRBl* 13:02 and D70 
were negatively associated with ACPA(+) RA in both SE-positive 
and -negative subjects (SE-negative: P= 0.0212, OR 0.59, 95% CI 
0.38-0.92; SE-positive: P= 0.0144, OR 0.53, 95% CI 0.32-0.87 
and SE-negative: P= 0.0011, OR 0.59, 95% CI 0.43-0.81: SE- 
positive, P= 0.0001, OR 0.56, 95% CI 0.42-0.75, respectively). 
D70 alleles other than *13:02 were protectively associated with 
ACPA(+) RA in SE-positive (P= 0.0076, OR 0.65, 95% CI 0.47- 
0.89), but not in SE-negative subjects. These data suggest that the 
negative association of D70 aUeles with ACPA(+) RA in SE- 
negative patients is mainly mediated by *13:02. Thus, the negative 
association of *13:02 with ACPA(+) RA was confirmed in SE- 
negative and -positive subjects. 

The DRB1*15:02 aUele was negatively associated with ACPA(-) 
RA in SE-negative people (P= 0.0008, OR 0.22, 95% CI 0.08- 
0.61). 167 and SI aUeles were negatively associated with 
ACPA(-) RA in SE-negative subjects (P= 0.0080, OR 0.45, 
95% CI 0.25-0.80 and P= 0.0008, OR 0.37, 95% CI 0.21- 
0.67, respectively). However, 167 alleles other than *15:02 or SI 
alleles other than *15:02 were not associated with ACPA(— ) 
RA. These data suggest that the negative associations of 167 and 
SI with ACPA(— ) RA in SE-negative subjects are mainly 
mediated by *15:02. Thus, the negative association of *15:02 
with ACPA(— ) RA was detected in SE-negative people. 



PLOS ONE I www.plosone.org 



6 



June 2014 | Volume 9 | Issue 6 | e99453 



DRBl^lSiOl Protects against RA 



< 

GC 



^ rN o^ m \D 



<— u-i u-i <— 



(N (N m ro 













o 


ro 




ro 


d 


^^^ 




ro 


o 


O 




VD 




o 


d 


d 


d 



rsi oj ro ro 

1 I I I 

LT) rs a* 



(-Ot/lUll/lLOrOLOLOU-lt/liyiLO 



LO LO LO l/l 



LO LO LO LO 



O t- <- 



ro vO <Ti 
r-. o rs o\ 
d r-^ d 



CTi v£) CO 



< 

0. 

u 
< 



O ^ O >- T- 



co rsi ro LO 



rM rsi rs 



uo vo 



-1.29) 


-29.04) 


-1.14) 






-0.86) 


-0.58) 


-1.15 


-1.37 


(0.84- 


(1.51- 


(0.40- 


(0.53- 


(0.47- 


(0.01- 


(0.31- 





o 


00 


d 




d 




O 


rs 


d 


d 


d 



4 4 4 





o 






00 




ro 














■=3- 






d 






d 


d 


ro 


LO 




a^ 


rN 






LO 






ro 




d 


d 


d 


ro 


d 


d 



m <— -— 



X O r- 

no ■— IT) 

o o 

d d 







LO 


1^ in 






rs 


rv. ,— 


CO 


LO 




>- 


CO 


m 


LO 




d 


ro 


d 


d 


d 


d 


d 


d rN 



LO LO LO LO 





00 


00 


o 


<- rM 




CO 


rs 








CO 






ro 






d 


d 


d 


d d 




d 


d 





LO 








d 










o 


CO 






o 




LO 




ro 


in 


d 


d 




d 


d 



+ 
< 

a. 

u 
< 



fvi 0^ m CTi 
vO ^ 
^ ro <— 



CO vO CO vO 



o 
u 























CO 






ro 


VD 














q 




O 
d 


ro 


CO 










ro 


in 




CO 




o 



(4.4) 


(2.8) 


ro 

d 


(5.6) 


LO 
OO 






in 





q 






OO 


vO 


ro 






ro 


in 


vO 


o 


CO 


LO 


in 
ro 


■=3- 
ro 
■=3- 


O 

in 



< 



u 
< 



to" ^ r-1 ^ ?M ?i R ^ 



O ^ 00 



^ ^ ^ 



s=r P 



























00 








(0.0) 


(3.6) 




(4.5) 


(0.9) 




(0.0) 


(3.6) 


(0.0) 


(5.5) 


(2.7) 


(0.0) 


13 (11 


15 (13 


(6.4) 


CO 


O 


■=3- 




Ln 






o 




O 




no 


O 







+ 
< 

u 
< 



_ ^ rj- „ 

^ i< fN "^J: 

vO O ' O 

^ LO fN ""^ 
^ «- 



^ O 

in 

o 



r-. — ^ 





CO 


CO 


G\ 




(N 


in 




in 


vO 


ro 






rsi 


in 





«- >- 



— — o ^ 
CO rN ^ ro 
>— ro ^ ^ 



rn 



CO CO CO 



~ ~ 



PLOS ONE I www.plosone.org 



7 



June 2014 | Volume 9 | Issue 6 | e99453 



DRB1*13:02 Protects against RA 



Then we examined the protective effects of *13:02 against 
ACPA(+) RA in the presence of predisposing alleles for ACPA(+) 
RA, DRB1*04:05 and *09:()1 (Table 5). As shown in Table 7, the 
risk for RA was decreased when these alleles were present together 
with *13:02 {*04:05: P= 0.0202, OR = 0.49, 95% CI 0.27-0.88; 
*a9;W: P=0.00,35, OR = 0.30, 95% CI 0.13-0.69). 

The protective effects of the * 15:02 allele against ACPA(-) RA 
were also analyzed in the presence of predisposing alleles (Table 7). 
DRB1*04:05 and *14:54 are potentially risk alleles for ACPA(-) 
RA (Table 5). The risk for ACPA(-) RA showed tendency 
towards decrease when these alleles were present together with 
'']5:02 {*04:05: P= 0.2665, OR = 0.39; *14:54: P= 0.6686, 
OR= 0.45), but these differences were not statistically significant. 

Certain amino acid residues in the HLA-DRp chain are 
associated with RA 

Finally, we analyzed the association with RA with respect to 
each amino acid residue in the HLA-DRP chain. Tyrosine at 
position 10 (lOY, P= 1.34x10"^", OR = 0.44, Pc = 4.59 x 10"'^ 
95% CI 0.37-0.52), serine at position 11 (US, P= 1.35x10"^", 
OR =0.44, Pc = 4.59xl0"'^ 95% CI 0.37-0.52), threonine at 
position 12 (12T, P= 1.35x10"^", OR = 0.44, Pc = 4.59 x 10"'^ 
95% CI 0.37-0.52), and aspartic acid at position 70 (70D, 
P= 1.15x10"^°, OR = 0.43, Pc = 3.91xl0"'-\ 95% CI 0.36- 
0.52) in the DRP chain showed strong protective associations with 
RA (Figure lA, open circles). Similar associations were observed 
with ACPA(+) RA (Figure IB), whereas aspartic acid at position 57 
(57D, P=0.0006, OR = 0.46, Pc = 0.0191, 95% CI 0.30-0.71) in 
the DRP chain showed a slight protective association with 
ACPA(— ) RA (Figure IC). Thus, association analysis suggested 
roles for specific amino acid residues in the HLA-DR(3 chain. 

Discussion 

Many groups have investigated associations between HLA- 
DRBl alleles and RA disease susceptibility. However, few studies 
have focused on protective effects of DRBl alleles against RA 
[1 1,12]. In the present study, we determined that the DRBl* 13:02 
allele plays a protective role in Japanese RA, especially in ACPA(+) 
RA, using RPE analysis (Table 2). A lower frequency of *13:02 
alleles in Asian patients with RA has been reported before 
[13,14,15,16,17]. In the genotype analysis, lower frequencies of 
the "HLA-DRB1*04:05/*13:02", or "*09:01 / *13:()2" genotypes 
in RA were observed (Table 4). Thus, the protective effects of 
*13:02 seem to overcome the predisposing effects of *04:05 or 
*09:01. Several studies have shown that certain DRBl alleles are 
negatively associated with RA and also some negatively associated 
allele groups defined by amino acid sequences, such as D70, 167, 
SI and DERAA (Table 3) [6,7,8,9,10]. Our results indicated that 
the protective effects of these allele groups were mainly 
attributable to *13:02 in Japanese RA, whereas they are 
attributable to * 13:01 in European RA [12]. The age at onset of 
*13:02 allele carriers was higher than non-carriers (Table S3), 
suggesting that the allele carriers of *13:02 may be associated with 
RA subsets with higher age at onset, and/or the age at onset may 
be delayed in the presence of the *13:02 allele. 

DRBl* 13:02 commonly belongs to the haplotype DRBl* 13:02- 
DQB1*06:04-DPB1*04:01, which shows evidence for positive 
selection in Japanese in recent history [27]. The DRB1*13:02 
allele is also a protective allele for cervical cancer [28], 
autoimmune hepatitis [29], and DPB1*04:01 is protective for 
hepatitis B infection [30]. Certain genes of this haplotype could be 
protective for these diseases, in addition to RA. 
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Table 7. HLA-DRBl genotype frequency in ACPA{+) and ACPA{-) RA patients and controls. 





ACPA{+) RA {n = 919) 


Control (n = 800) 


P 


OR 


95%CI 


*04:05/'13:02 


28 (5.9) 


21 (11.4) 


0.0202 


0.49 


(0.27-0.88) 


'09:01/*! 3:02 


8 (3.2) 


21 (9.9) 


0.0035 


0.30 


(0.13-0.69) 


*04:0U* 13:02 


3 (4.6) 


2 (11.8) 


0.2755 


0.36 


(0.06-2.37) 


*01:0U* 13:02 


4 (2.7) 


3 (3.6) 


0.7061 


0.75 


(0.16-3.44) 


* 14:54/* 13:02 


0 (0.0) 


6 (13.3) 


0.0263 


0.07 


(0.00-1 .28) 




ACPA(-) RA (n = 110) 


Control (n = 800) 


P 


OR 


95%CI 


*04:05/* 15:02 


2 (5.3) 


23 (12.4) 


0.2665 


0.39 


(0.09-1.74) 


*09:01/'15:02 


1 (3.3) 


28 (13.1) 


0.2228 


0.23 


(0.03-1 .74) 


*04:01/* 15:02 


1 (25.0) 


1 (5.9) 


0.3524 


5.33 


(0.26-110.80) 


*01:01/*15:02 


0 (0.0) 


5 (6.0) 


1.0000 


0.68 


(0.04-13.19) 


'14:54/*15:02 


1 (7.7) 


7 (15.6) 


0.6686 


0.45 


(0.05-4.06) 



ACPA: anti-citrullinated peptide antibody, ACPA(+): ACPA-positive, ACPA(— ): ACPA-negative, RA: rheumatoid arthritis, SE: Shared epitope, OR: odds ratio, CI: confidence 
Interval, Allele carrier frequencies are shown in parenthesis (%). Association was tested by Fisher's exact test using 2 x2 contingency tables. Upper row: "B/*13:02" vs. "B/ 
other than '13:02". Lower row: "Bl'15:02" vs. "fi/other than '15:02". 
doi:l 0.1 371/journal.pone.0099453.t007 



It was reported that SE alleles are strongly associated with 
ACPA(+) RA, but weaHy with ACPA(-) RA [1], and this was 
confirmed in the present study. We documented protective effects 
of DRB1*13:02 against ACPA(+) RA and DRB1*15:02 against 
ACPA(-) RA in Japanese. Although the sample size of ACPA(— ) 
RA is not large enough, the protective effect of *15:02 against 
ACPA(— ) RA was also reported in another study [31], supporting 
the results. These findings could be explained by differences in the 
pathogenesis of ACPA(+) and ACPA(-) RA. Although the 
genotype of DRB1*03/ *13 was reported to be associated with 
ACPA(— ) RA in a European population [1 1], such an association 
was not found in the current study. 

Amino acid residues lOY, 1 IS, 12T, and 70D of the HLA-DRP 
chain were negatively associated with RA (Figwe 1 A). Amino acid 
residues 1 1 and 70 form the HLA-DR peptide-binding groove 
[32]. These data suggest the involvement of peptide antigens 
bound to specific HLA molecules in controlling the development 
of RA. Associations of amino acid residues 10, 11, 12, 13, 33, 37, 
47, 67, 70, 96 and 98 of the HLA-DRP chain were reported in 
European ACPA(+) RA [33], showing slightly different association 
pattern from the results of this study (Figure IB). However, 
associated amino acid residues 10, 11, 12, 13, 33, 57, 70, 96 and 



98 of HLA-DRP chain in Korean ACPA(+) RA [33] were more 
similar to the results (Figure IB), reflecting the difference of DRBl 
allele frequencies between European and Asian populations. 

The negative association with the DRBl * 13:02 allele needs to be 
confirmed in future independent studies. Because the distribution 
of HLA alleles in other ethnic populations is different from the 
Japanese, the protective role of some DRBl alleles in RA in other 
populations should be determined. 

Thus, the present study identified a negative association of 
DRBl* 13:02 with Japanese RA; our findings support the 
protective role of DRBl* 13:02 alleles in the pathogenesis of 
ACPA(+) RA. 

Supporting Information 

Table SI HLA-DRBl homozygous frequency in the RA 
patients and controls. 

(PDF) 

Table S2 HLA-DRBl allele frequency in the RA patients 
and controls. 

(PDF) 
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Figure 1. Associations of amino acid residues in the DRp chain with RA (A), ACPA-positive [ACPA(+)] RA (B), and ACPA-negative 
[ACPA(— )] RA (C). Corrected P (Pc) values were calculated by nnultiplying the P value by the number of amino acid residues tested. Associations 
were established by Fisher's exact test using 2x2 contingency tables. Positive associations are indicated by filled circles and negative associations by 
open circles. 

doi:1 0.1 371/journal.pone.0099453.g001 
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Table S3 Age at onset of HLA-DRBl allele carrier or 
non-carrier in the RA patients. 
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